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In order to successlully describe DIS data, one must take heavy quark mass ef- 
fects into account. This is often achieved using so called variable flavour number 
schemes, in which a parton distribution for the heavy quark species is defined 
above a suitable matching scale. At small x, one must also potentially include 
high energy corrections to this framework arising from the BFKL equation. We 
outline the definition of a variable flavour scheme which allows such corrections to 
be consistently implemented alongside a fixed order QCD expansion. Results of 
a global fit using this scheme are presented. We also discuss an extension of the 
resummation to NLL order. 



1. Variable Flavour Number Schemes 

In DIS one can produce heavy quarks in the final state via boson gluon 
fusion. Taking into account all the relevant Feynman diagrams, one may 
absorb the mass effects into the coefficient functions in a chosen factori- 
sation scheme. This is a fixed flavour number scheme, and is undesirable 
at high Q 2 in that the coefficients diverge due to collinear singularities. 
Instead, one may define a heavy quark distribution above a suitable match- 
ing scale (here chosen as Q 2 — M 2 ), which evolves according to DGLAP 
equations. This is a so-called variable flavour description, as u/ changes 
across Q 2 = M 2 . Equivalence of the FF and VF descriptions imposes a 
consistency relation between the partons above and below the matching 
scale 1 : 
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which defines the heavy matrix elements {Aij}. This leads to an ambiguity 
in the VF coefficient functions. For example, at 0(as), one finds: 

cZ$ 0) ®aV + c%v = c[JV (2) 

Provided this relation is satisfied, one has a valid VF number scheme. How- 
ever, one is free to shift terms <~ 0(M 2 /Q 2 ) between the VF coefficients. 
Each choice corresponds to a different VF number scheme, and this trans- 
formation does not change the collinear factorisation scheme. In practice 
one specifies the placement of M 2 /Q 2 terms by fixing one of the VF coef- 
ficients according to some choice. 

2. Heavy Quark Production at Small x 

Our aim is find a VF number scheme that has the same definition at fixed 
order as it does in the high energy expansion (small x). It is convenient 
to adopt the DIS factorisation scheme for this purpose, as then the LL 
impact factors coupling the virtual photon to the BFKL gluon ladder via a 
massless quark pair are easily interpreted in terms of the splitting function 
P qg and longitudinal coefficient Cl 9 - There is no mixing between Ci g and 
P qg as the former quantity is zero in this scheme. One must also choose the 
placement of M 2 /Q 2 terms. Previously, the ACOT group has proposed 
the identification C^^ 0>>VF = 5(1 — x') (the ACOT(x) scheme 2 ), where 
x' is the scaled Bjorken variable x' = x(l + AM 2 /Q 2 ) taking into account 
the kinematic constraint at the heavy quark vertex. We propose the same 
identification in the DIS scheme n.b. C 2 h = 5(1 — x') to all orders. This 
then specifies completely C%g and C% F and we call this the DIS(x) scheme 
by analogy with ACOT. It is then easy to interpret the heavy quark impact 
factors at small x, and in double Mellin space one finds the resummed 
quantities: 

A Hg ( 7 , N, Q 2 /M 2 ) = h 2 ( 7 , N, Q 2 /M 2 )\ Q2 ; (3) 

M 2 

C\ F g = h 2 (j, N, Q 2 /M 2 ) - A Hg {l, N, Q 2 /M 2 ); (4) 
Cl F g =C[ F g =h L ( 7 ,N,Q 2 /M 2 ), (5) 

with N Mellin conjugate to x and 7 conjugate to Q 2 /A 2 . Results in x and 
Q 2 space depend upon how one solves the BFKL equation. 

3. A LL Resummed Global Fit with Running Coupling 

We implemented the DIS(x) scheme in a LO global fit with LL 
resummations 3 , where the solution of the BFKL equation also included 
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the running coupling 4 . Although this latter effect is strictly subleading, 
the resulting resummcd coefficients show a considerably softer divergence 
at small x. The resummed prediction is compared with data alongside a 
standard NLO prediction from a global fit in figure f. The resummcd 
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Figure 1. Comparison of a LL resummcd fit (with running coupling) to data from 
HERA and NMC. Also shown is a NLO fit result. 



prediction is excellent at small x, whereas the NLO description starts to 
fail because its slope is insufficient as Q 2 increases. However, a problem 
with the resummed fit (not shown) is that the predicted gluon and sea dis- 
tributions (hence Fl and jet cross-sections) arc too low in the moderate 
x region. This indicates a need for NLO as well as NLL corrections, to 
suppress resummation effects to lower x values. 

4. Toward a NLL Analysis with Exact Gluon Kinematics 

In calculating the impact factors coupling the virtual photon to the BFKL 
gluon ladder, one assumes that the longitudinal momentum fraction of the 
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gluon at the top of the ladder is the same as Bjorken x (corresponding to the 
bottom of the ladder). This is not true beyond LL order, and the impact 
factors with the exact gluon kinematics have been previously calculated . 
We have shown that they seem to contain most of the NLL correction to the 
impact factor 6 , and thus pave the way for an approximate NLL analysis of 
scattering data. A preliminary result for the NLL corrected P qg is shown 
in figure 2. The NLL corrections soften the low x divergence below the LL 
prediction and delay the onset of asymptotic behaviour to lower x values. 
This is hopefully what is needed to obtain a good description of scattering 
data over the whole of the accessible x range! 
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Figure 2. Preliminary result for P qg with NLL and running coupling corrections. Also 
shown are LL (with running coupling) and fixed order results. 
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